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Abstract 

The problem of separation of variables in a dyon-dyon system is discussed. A 
linear transformation is obtained between fundamental bases of this system. Com- 
parison of the dyon-dyon system with a 4D isotropic oscillator is carried out. 
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1 Introduction 



In this paper, we have calculated the matrix between the spherical and parabolic bases of 
a dyon-dyon system [1] belonging to the same energy level. This matrix is a generalization 
of the Park-Tarter matrix known from the theory of hydrogen atom [2, 3] to the case when 
the Coulomb center carries not only the electric but also magnetic charge. Like the Park- 
Tarter matrix, our matrix is expressed through the Clebsch-Gordan coefficients C^2^;6/3' 
however, in our case a ^ b, in contrast to the case of a hydrogen atom. We have also traced 
the connection of the dyon-dyon problem with that of a 4-dimensional isotropic oscillator. 
As is known [4], these problems are related to each other by the Kustaanheimo-Sticfcl 
transformation [5] supplemented with the 4th (angular) coordinate. We have shown that 
the coefficients C^a;6/3 coincide with the ones [6] of the expansion of the double polar basis 
over the Euler basis of a 4-dimensional isotropic oscillator. 

2 Dyon— Dyon System 

A dyon-dyon system in the space IR^ is described by the equation 



dx-i he ^ I 



^ + ^h^ + -)^ = (1) 




where 



r(H — Xg) 

and s — eg/hc — 0, ±1/2, ±1, .... Each value of s describes its particular dyon-dyon 
system. At s = 0, eq.(l) is reduced to the Schrodinger equation for a hydrogen atom. 
When s 7^ 0, equation (1) preserves 0(4)-symmetry and therefore variables in it are 
separated into spherical, parabolic, and prolate spheroidal coordinates [1]. 

The system (1) possesses a singularity on the axis 2:3. It is also possible to consider 
systems with singularities either on the semiaxis 2:3 > or on X3 < 0, i.e. they are 
described by the vector potentials 

Af= , ^ {tx2,±x,,Q) 
■> r{r^x^) 

and are connected with the system (1) by the gauge transformations 

Af^ = A, + V'(±)(f) = ^(f) exp 

with the gauge function f^^^ — ±2g aictcLnx2/ Xi. 

The variables in eq. (1) are separated in spherical and parabolic coordinates. 
In the spherical coordinates 

Xi — r sm.6 cosip, X2 = r sn\6 snap, X3 = rcos9 (2) 
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the wave function of the dyon-dyon system is of the form [7] 



where the functions zjfj^{6) and Rnlmi''^) normahzed by the condition 



is) 



smezit{e)zi:^{e)de = 6,,k, 



dr = 1 



are given by the formulae 

ztliO) = Niiil - cos^)^(l + cos^)^Pil— l>-l)(cos^) 



, . \m-3\ + \m+3\ 

/ 2r 2 



ronj \ronJ 



„ I , \m — s\ + \m + s\ , , , , , 2r \ 
F\-n + k+ ^—^ + 1; 2A; + m - s + m + s + 2; 

Here P^"'^) (x) are Jacobi polynomials; tq — fi' jM^e} is the Bohr radius. The normaliza- 
tion constants A^^^^ and C^^^ equal 



(2A; + 


m — s 


+ |m + s| + + 


m — i 


;| + |m + s )! 


2|m-. 


5|+|m+s|+ir(A; + 


m — s 


+ l)r(A;+ 1? 


Ti + S] + 1) 



1/2 



(^) 

nkm 



1 



(n + A; + 



|m— s| + |m+s 
2 



7^2^3/2 + |m - s| + |m + s| + 1)! (n - k - 1^-^1+1"^+^! _ i 

Quantum numbers run over the values n — 1, 3/2, 2, A; = 0, 1, ■■■kmax, where 

\m — s\ + \m + s' 



The energy spectrum of the system is of the form 
In the parabolic coordinates 

Xi = 

upon the substitution 



- 1 



^7]COS(p, X2 = \/^7]Sm(p, X3 = -{^-7]) 



(3) 



(4) 



27r 



the variables in (1) are separated, which results in the system of equations 

/i = 



di v dd 



- (v—) + 

dr] \ dr] I 



2h 



2 S 77 r Pi 
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2 V 



where 



4e 
477 

Moe^ 



+ /32 



/2=0 



A + /92 = ^ 



(5) 



At s = 0, these equations coincide with the equations for a hydrogen atom in the parabolic 
coordinates [8], and consequently, 



3/2 Jni,m+s\'i ) Jn2,m-s\V ) 



where 



lT{p+\q\ + l) 







exp 



27r 



r(kl + i)V p! 

Here rii and 712 are non-negative integers 



Ir^nj \rQnJ \ Ton J 



X \ 



ni 



|m + s| + 1 



+ 



h 



V-2Moe* 



|m - s| + 1 

Pi, ^2 = z h 



h 



/32 



from which and (3), (5) it follows that the parabolic quantum numbers ni,n2,m and s 
are connected with the principal quantum number n as follows: 



\m — s\ + \m + s 
n = ni + ^2 H h 1 



(6) 



3 Park— Tarter Generalized Matrix 

We write the searched expansion in the form 



(7) 



k=0 



Our purpose is to calculate the coefficients T^f^jfcm; i-^- ^^e Park-Tarter generalized 
matrix. The usual Park-Tarter matrix is the matrix T^^l^j^^ at s = 0. 
We substitute 

^ = r(l + cos 6*), T] = r(l — cos^), 
into the left-hand side of expansion (7), let r tend to infinity, take the formula 

F(-n; c; x) ~ (-l)"f^^^^ - 00) 



3 



and the orthogonality condition for the function zj^J^ into account. All this leads to the 
formula 



rpi^) _ (_-\)kiD(s) t{s) 



where 



B, 



is) 

ninzkm 



{2k + 


m — s 


+ m + s 


+ l)k\{k + 


m — 6 
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m + s )! 


22n+\m-s\+\m+s\Y{k + 


m — s 


+ l)r(A;+| 


m + + 1) 



'[n-k- I'^-^l+l'^+^l -i)! (n + A;+ ["'-^l+l'^+^l )! 
(ni)!(n2)!r(ni + |m + s| + l)T{n2 + |m — s| + 1) 

and the second factor is equal to the integral 

-1 

Then taking advantage of the Rodrigues formula [9] 



1/2 



\fS+n 



and the integral representation for the Clebsch-Gordan coefficients [10] 

(2c+l)(J+l)!(J-2c)!(c + 7)! 



(J - 2a)!(J - 2b)\{a - a)\{a + a)l{b - p)\{b + /3)\{c - 7)! 



1/2 



_l-^a-c+/3 



/(I - x)— (1 + xf-^-^ [(1 - xy-"^{i + x) 



2J+1 



\J-2b 



dx 



{J — a + b + c), we obtain 



nk 



where 



ni + n2 + \m + s\ 



rii + n2 + \m — s\ 



6/3 



, c = k + 



\m — s\ + \m + s\ 



a — 



rii — n2 + \m + s\ 



n2 — rii + \m — s\ 
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\m — s\ + \m + s\ 



2 ' ' 2 ' ' 2 

At s = formula (8) turns into the Park-Tarter formula, as would be expected. 
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4 Dyon— Dyon System and 4D Oscillator 

Let us demonstrate that if in eq. (1) we make the changes 



d 



(9) 



(7 e [0, 47r)), it will transform into the Schroedinger equation for a 4D isotropic oscillator. 
Equation (1) in the spherical coordinates is of the form 



1 A (r-^] + 1 
dr \ dr J 

2is cos 9 dijj 



sin^ 6 dip sin^ 6^ ' fi 
Prom (9) and (10) we have 



1 d f . djj\ 1 ay 

'^de\ '89 I ^ sin' e dip'' 



, 2Mo / , en , ^ 



(10) 



1^ r'^ 
r' 9r \ 9r 



^2 



2Mo 



where 



1 d f . 1 



sin /? dp 



dpi sin' p\da^ 



d' . 92 
2 cos /^T^^ + 



dad^ ^ 



Here we change the notation: P — 9 and a — (p. If we now pass from the coordinates 
r, a, P, 7 to the coordinates 



P -42+2. . . P 

Uq + lUi — u COS — e 2 ^ + lu^ — u sm — e 2 



(12) 



with m' = r, take into account that 



^^^d_(,d_^_A_j, 
du^ du \ du ^ 



and introduce the notation 



E = 4e^ 



then equation (11) will turn into the Schrodinger equation for a 4D isotropic oscillator 

^(m) = 



9' ^2Mo('^ Mou;2^2> 



^ V 2 ^ 

whose energy spectrum is given by the formula 

En = huj{N + 2) 



(13) 
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Introducing the double polar coordinates 

uo + iui ^ Pie^"^\ U2 + ius ^ P2e"^^ (] 
from formulae (2), (4), (12), and (14) we get the relations 

^ = 2pl, T) = 2pl, (P = (P1 + (P2, 1 = ^1-^2 

which lead to the formulae 

i>NJMiM2{u, a, 7) = 'insj^5„ 6j^j im^-m^Wm^+m^I Sm.U^ Ss.M^i^nkmsjr, 0, Ifi, 7) 
i'NiN2mim2{Pl^ P2, ^1, ^2) = 4nW — . Afi (^no .Afo m^+^m^ (5„^ -m^ jln, nom.'^ '^i l) 

generalizing the earlier results [6, 11]. 

Now we are able to write the expansion [6] 

N/2 

V'iViiVamima (Pl, P2, ^1, ^2) = Yl ^mN^^m^'^NJM.M^ {% a, /3, 7) (] 



J — Jfr 



where 



TTrNJMiM2 _ i7r<I>^co,70 (-1 
'^^^ NiN2mim2 ^ "-^ao,ao;6oA ^ 

N + Iniil — \m2\ , A^— |mi| + |m2| 
oo = ^ , bo = , Co = J 

N + \mi\ - \m2\ „ A^-|mi| + |m2| + |m2| 

tto = z -'V2, Po = iVi, 7o 



4 ^' 4 " 2 

The lower limit of summation in (15) and quantity $ are given by the expressions 

Jmin = ^ i\Mi - M2I + |Mi + M2I) 

^ _ Ar , T l"^il + l"^2| m2 + \m2\ 
^ 2 2 

We conclude with the following two comments: 
(a) Using formulae (2) and (12) and considering that r — u^,9 — P,(p — a, one can eas 
show that 

Xi = 2{uoU2 + UiUs) 
X2 — 2(moM3 — U1U2) 

— Uq + ul — ul — ul 



i. {uo + iui){u2 + ius) 

7 = — in — 

2 {uq - iui){u2 - ius) 
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The first three fines are tfie transformation ]R"^ ]R'^ suggested by Kustaanfieimo and 
Stiefel for the regularization of equations of celestial mechanics [5]. Later, this trans- 
formation found other applications, as well [12, 13]. This transformation supplemented 
with the coordinate 7 was used for the "synthesis" of the dyon-dyon system from the 4D 
isotropic oscillator [4]. 

(b) It is known [6] that diagonal (mi = 7712) elements of the matrix W^;'^^l,^%,^ with N 
even coincide with the Park-Tarter matrix. From formula (16) it follows that the remain- 
ing elements of the matrix W^^^^^J^^^^ have also a physical meaning: these are elements 
of the generalized Park-Tarter matrix for the dyon-dyon system. 



5 Degeneracy of the Energy Levels 

Let us discuss the problem of multiplicity of degeneration of tfie energy levels (3) and 
(13). From formula (6) it follows tfiat at fixed n, m and s tfie energy levels are degenerate 
witfi tfie multiplicity 

, |m — s| -|- |m -|- s| 



For s > tfie multiplicity of degeneration of levels (3) at fixed s and n is 

\m\>[s] \m\<[s]-l 

wfiere tfie upper limit of summation is determined from tfie condition g^^ > 0, 

|m — s| -|- |m -|- s| < 2n — 2 

Tfierefore, 

9n- E {n-s) + 2 J2in-m)^{n-s){n + s) (17) 

m=— [s]+l m=[s] 

Tfie same result follows from analogous computations also wfien s < 0. 

Tfie quantum numbers s and n in formula (17) assume simultaneously eitfier integer 
or fialf-integer values, and tfius, we fiave 

9n^ E ^^ = ^n(2n-l)(2n + l) 

s=-n+l 

wfiere gn stands for tfie multiplicity of degeneration of tfie energy levels (13) of tfie 4d 
oscillator. Since N — 2n — 2, we arrive at tfie known result 

gN^l{N+l){N + 2){N + 3) 
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